High-order optical nonlinearity at low light levels 



O 

(N 

O 

Q 

o 



o 
q 

O 
>V 

^ : 

Oh. 



(N 

> 

O 

o 



X 



Joel A. Greenbergj and Daniel J. Gauthier 

Department of Physics and the Fitzpatrick Institute for Photonics, Duke University, Durham, NC 27708, USA 

(Dated: January 4, 2012) 

We observe a nonlinear optical process in a gas of cold atoms that simultaneously displays the 
largest reported fifth-order nonlinear susceptibility x' 5 ' = 1.9 x 10 -12 (m/V) 4 and high transparency. 
The nonlinearity results from the simultaneous cooling and crystallization of the gas, and gives rise 
to efficient Bragg scattering in the form of six- wave-mixing at low-light-levels. For large atom-photon 
coupling strengths, the back-action of the scattered fields influences the light-matter dynamics. This 
system may have important applications in many-body physics, quantum information processing, 
and multidimensional soliton formation. 
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Since the first observation of second-harmonic gener- 
ation of a ruby laser beam over 50 years ago, there has 
been sustained effort to improve the nonlinear optical 
(NLO) interaction strength of materials. One ultimate 
goal is to realize nonlinear interactions at the single- 
photon level, which will lower the operating power of 
devices. There is also a growing need for single-photon 
nonlinearities for quantum information applications 

The strength of a material's nonlinear optical response 
is characterized by the n-th order susceptibility tensor 
> which relates the nonlinear polarization of the ma- 
terial P NL to the electric field strength E via P NL = 
e [V (3) : EEE + V (5) = EEEEE + ...] for an isotropic 
material 0] ■ Because the magnitude of the susceptibility 
typically decreases with increasing order, most low-light- 
level studies to date have focused on lower-order pro- 
cesses. For example, there have been numerous observa- 
tions of low-light-level NLO interactions based on third- 
order (x^) processes created via electromagnetically- 
induced transparency (EIT) where a strong cou- 

pling beam creates a quantum interference effect that 
simultaneously renders the medium nearly transparent 
while enhancing the nonlinearity Q • Other recent obser- 
vations of strong third-order NLO effects include a two- 
photon absorptive switch actuated using < 20 photons 
interacting with atoms in a hollow-core photonic bandgap 
fiber [sl and an optical pattern-based switch using ^600 
photons interacting with a warm atomic vapor [9j . 

Despite the success of these approaches, some applica- 
tions require or can benefit from higher-order nonlinear- 
ities. Materials with a large, fifth-order (x ) response, 
for example, can lead to new phenomena, such as liquid 
light condensates and transverse pattern and soliton 
formation and play an important role in quantum 
information networks through enabling 3-qubit quantum 
processing 112 . 1 1 31 ] , providing new sources of correlated 
pulse pairs |l4j |. and acting as quantum memories (lilllij . 
Quintic media can also improve high-precision measure- 
ments and reduce phase noise for enhanced interfer- 
ometry performance [171 ]. Thus, the realization of effi- 



cient materials is important for fundamental studies 
in quantum nonlinear optics as well as improving the per- 
formance of NLO devices. 

In this Letter, we report the discovery of a dissipation- 
enhanced NLO process that gives rise to the largest fifth- 
order (x ) NLO susceptibility ever reported while si- 
multaneously having high transparency. Our NLO mate- 
rial consists of a gas of cold atoms initially in thermal 
equilibrium, which is illuminated by weak, frequency- 
degenerate laser beams (frequency lu). The light fields 
(both applied and self-generated via wave mixing) act on 
the atomic center-of-mass motion to cool and crystallize 
the gas and lead to Bragg scattering via the generated 
atomic density gratings. This cooling arises via the dis- 
sipative Sisyphus force [3] and occurs efficiently (i.e., 
requires the scattering of only tens of photons per atom) 
so that absorption can be made small. Thus, in contrast 
to previous studies of wave mixing via atomic bunching 
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that produce a third-order NLO response 
sipative effects are crucial to the nonlinearity discussed 
here and cause the lowest-order nonlinearity to be fifth- 
order in the applied fields. 

Surprisingly, the achievable nonlinear interaction 
strength observed in our experiments from the x^ 5 "* re- 
sponse is just as large as that obtained in previous ex- 
periments dominated by a x response. This strong 
light-matter coupling enables the scattered fields to al- 
ter the atomic spatial organization and momentum dis- 
tribution, resulting in greatly reduced group velocities, 
enhanced atomic coherence times, and optical instabil- 
ities. Because of this back-action, our system may be 
particularly interesting for studies of strongly-correlated 
many-body physics with long-range interactions (22) . 

To make our discussion concrete, we consider the situ- 
ation shown in Fig. [I] where optical fields interact with 
a cloud of cold atoms in a pencil-shaped geometry. A 
pair of strong, counterpropagating pump fields (intensity 
I p , wave vectors ±k p ) are inclined by an angle 0=10° rel- 
ative to the cloud's long axis and are linearly-polarized 
with orthogonal polarizations (linJ_lin configuration). A 
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weak signal field (intensity / s , wavevector k s ) is injected 
along the z-axis with a polarization parallel to the nearly- 
countcrpropagating pump field. An idler field (Jj, -k s ), 
self-generated via NLO wave mixing, counterpropagates 
with and has a polarization that is orthogonal to the sig- 
nal field. We refer to the signal and idler fields as probe 
fields. 
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FIG. 1. Schematic of the experimental setup. We take z — 
(z — L) to be the left (right) side of the cloud. 

In our experiment, we use an anisotropic magneto- 
optical trap (MOT) to confine 87 Rb atoms in the 
5 2 Si/ 2 (F = 2) state within a cylindrical region of length 
L=3 cm (along z) and diameter 14^=300 fxm. (along x 
and y) [23J. This configuration enables us to achieve 
typical atomic densities of ~ 5 x f0 10 cm 
isotropically cooled to T e9 =20-30 /jK. The 
probe beams are tuned below the 5 2 Sx/2{F = 2) — > 
5 2 P3/2(F' = 2) transition (transition frequency W23) by 
|A|/r = 3 — 20 (A = u — W23) and have diameters of 
3 mm and 200 /im, respectively. We use pump beam 
intensities of up to a few mW/cm 2 and a signal beam 
intensity of 3 /zW/cm 2 . Thus, all beams are well below 
the off-resonance electronic saturation intensity i~ t = 
I sat [l + (2A/r) 2 ], where I sat = Ae ch 2 / (fi 2 T 2 ) = 1.6 
mW/cm 2 is the resonant saturation intensity, T/2tt = 6 
MHz is the natural transition linewidth, /1 = 2.53 x 10~ 29 
C-m is the reduced transition dipole moment, eo is the 
permittivity of free space, and c is the speed of light in 
vacuum. 

To measure the system's NLO response, we cycle be- 
tween a wave-mixing and a cooling and trapping phase. 
During the wave-mixing phase, we turn on only the 
pump and signal beams for ~ 1 ms and record the time- 
dependent intensities of the signal [I s (z — L)] and idler 
[h(z = 0)] beams as they exit the cloud. After ~ 200 /us, 
the system reaches a steady state that persists for ~lms 
(at which time expansion of the cloud in the y-direction 
reduces the density and, consequently, the nonlinearity). 
We then cool and trap the atoms with only the MOT 
beams for 99 ms and repeat the cycle. We have verified 
that the MOT magnetic fields, which remain on during 
the experiment, do not affect the NLO response 
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Figure [2^i shows the steady-state reflectivity R = 
7j(0)/7 s (0) as a function of I p for different A, where 
the points (solid curves) correspond to experimental 
measurements (theoretical predictions, which we discuss 
later). The reflectivity scales superlinearly with I p and 



I sa t)- Beyond R ~ 1, an optical instability occurs that 
gives rise to self-generated signal and idler fields in the 



absence of an injected signal beam [24J. This instability, 
coupled with the rapid variation of R on I p , ultimately 
limits the largest values of R that we can measure. 
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FIG. 2. a) Dependence of R on I p for A/r = -3, -5, -7.3, -12.7, 
and -18.8 (from left to right, respectively), b) Signal-beam 
transmission for I p = 0. c) Nonlinear phase shift correspond- 
ing to the data in a) obtained via Eq. 4. The vertical dashed 
line indicates the value of the characteristic intensity I c . The 
inset shows that, for I p <C Ic, ft is well-fit (r 2 =0.994) by a 
quadratic function (where we show results for A/r = —3). 
d) Cross-phase modulation figure of merit as a function of 
I p for the same detunings as in a). The fact that all points 
follow the same curve indicates that both j3 and «a scale as 
(A/r) -2 . For all figures, oa=oL = 13.4(±0.5) (uncertainty 
corresponds to one standard deviation), the points correspond 
to the experimental data, and the solid curves correspond to 
the theory via Eqs. [3] and [4] 

For a fixed I p , R decreases with increasing |A|. Nev- 
ertheless, this decrease occurs sufficiently slowly that we 
can still obtain large R while operating far from reso- 
nance (where absorption is minimal). Figure shows 
the signal-beam transmission T = I S (L) / I s (0) in the ab- 
sence of the pump beams, given by T = exp(— a&.L), 
where a& = rjhhjT /2I^ at is the off- resonance absorption 
coefficient and r\ is the average atomic density. For all 
of the detunings shown in Fig. [5^, T > 0.65 and ap- 
proaches 1 for the larger detunings. Thus, this light- 
matter interaction allows us to realize large nonlineari- 
ties with high transparency for low input intensities (e.g., 
we measure R = 1 for I p = 1.5 mW/cm 2 by working at 
A = -5r where T = 0.87). 

We interpret the mechanism underlying the probe 
beam amplification as Bragg scattering of pump photons 
into the probe beam direction via an atomic density grat- 
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ing. The forces resulting from the interference of a pump 
and nearly counterpropagating probe beam generate an 
atomic density modulation along G — k p + k s that is 
phase-matched for Bragg scattering. The resulting non- 
linear atomic polarization is P NL = rnj, 2 bE p /hA, where 
b is the amplitude of the density modulation along G and 
E p is the pump electric field [25| . Thermal motion limits 
the degree to which atoms can be localized along G so 
that lower temperatures result in larger values of b and, 
consequently, higher Bragg scattering efficiencies. 

Because of this connection between atomic temper- 
ature and NLO scattering efficiency, we find that the 
dissipative optical lattice formed by the pump beams 
plays an important role in the nonlinearity (in contrast 
to other works in which the pump-beam lattice does not 
directly enhance the light- matter coupling [26|, The 
pump-beam lattice causes Sisyphus cooling, which trans- 
forms the atomic momentum distribution along k p from 
a Maxwell-Boltzmann distribution to one that is well- 
described by a double-Gaussian function [28l . [29j j . One 
interprets the gas as a non-thermal system consisting of 
a cold component of atoms well-localized in the pump- 
pump lattice at a temperature T' c and an unbound, hot 
component at temperature T' h undergoing anomalous dif- 
fusion. The cooling process transfers atoms from the hot 
to the cold component, where the fraction of atoms in 
each component is fh.c, respectively. In our experiment, 
we find that T' c = 2.5 fj,K and T' h ^ 25 /xK and that both 
are largely insensitive to I p and A [241 ] . 

The pump-beam lattice assists in cooling and loading 
atoms into the lattice along G, which is phase-matched 
for scattering pump light into the probe beam directions. 
For weak bunching and (2A/r) 2 > 1, we find that [Io| 



hT e c(E s E* p + E*E p ) / /, 



4(A/r) 



knT r 



fh 
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(1) 



where E s (Ei) is the signal (idler) electric field strength 
and T c ^h are the temperatures along G (which are ap- 
proximately equal to T' c h for 9/2 <C 1). We solve for f c ^ 
by calculating numerically the steady-state momentum 



distribution for a J„ 



-> J e = 3/2 transition using 



of I c on A plays an important role in achieving large NLO 
coupling strengths in our system. 

Substituting the atomic polarization given above into 
Maxwell's equations, we find that the steady-state cou- 
pled wave equations for the signal and idler beams be- 
come 



dE. 
dz 



-t^ = iKE s + ifiE*, 



dEl 
dz 



in* E* 



if3E s 



(2) 



where we make the constant-pump-beam approximation 
and assume that the optical fields adiabatically follow 
the evolution of the atomic density grating. We define 
the nonlinear coupling coefficient 



P = P c + P h 



fc 



fh 



«A=o^r i p 

64(A/T) 2 1Z \k B T c ' k B T h 



(3) 



such that j3L corresponds to the nonlinear phase shift 
imposed on the weak probe beams by the interaction and 
k = n R + iKi = P + a A =o(-2A/r + i)/8(A/T) 2 . 
From Eq. [2j we find that [33| 



R 



|^sin(wL)| 5 



\wcos(wL) + K/sin(u>-L)| ' 



(4) 



where w — (\P\ 2 — k 2 ) 1 / 2 . The solid lines in Fig. [2] 
demonstrate that the predictions of Eq. [4] fit the exper- 
imental data well (Xred = 1-6) over all measured values 
of A and I p . Using the measured values of aA=oI, A 
and T c ^, the only free parameter used to obtain this fit 
is I c , which we assume is independent of A. We find that 
I c = 1.5(±0.5) mW/cm 2 . We can therefore access NLO 
interaction strengths comparable to those of resonantly- 
driven atoms while the atoms are nearly transparent, 
since I c = I sa t even for |A/r| 3> 1. 

In order to quantify the nonlinear response, we use Eq. 
[4] to determine P directly (see Fig. [2}:). For small pump 
intensities (i.e., I p <C I c ), a- Taylor series expansion of p 
about I p = yields 



a ^ "A^o^r I p ( I p /I c 1 - I p /I c 
P ~ 64(A/r) 2 / sat \k B T c k B T h 



(5) 



a Bloch state approach 1^ and find that f c initially in- The quadratic (linear) components of p correspond 



creases linearly with I p before asymptoting to 1. While 
no analytic solution for f c exists, the simple functional 
form f c = l—ff l = tanh( Jp/J c ) fits the numerical results 
well, where the characteristic intensity for the cooling 
process (and, therefore, the nonlinearity) is I c . We find 
that I c is directly proportional to and of the same or- 
der as the so-called decrochage intensity Id, which is the 
intensity where the root mean squared width of the mo- 
mentum distribution is a minimum ( i.e., where cooling 
optimally balances diffusive heating) [18(. Perhaps sur- 
prisingly, both the cooling and heating rates vary with 
detuning in such a way that Id and, therefore, I c is inde- 
pendent of A [32|. As we discuss later, the independence 



to a x 



(■->) 



(x (3)) 



according 



to 



response 

\k p \/ (4Leoc){xx%yIp + (2e c)~ 1 Xxxxxyylp] for the polar- 
izations shown in Fig. [T] Using this relationship, we find 
XxLxyy = 1.9(±0.3) xlCT 12 (m/V) 4 for A = -3I 1 (see in- 
set in Fig. [5b). This x^ response is due almost entirely 
to the cold component (since T c /Th = 0.1) and gives 
rise to six- wave mixing (SWM), where the pump beams 
simultaneously undergo Bragg scattering and transfer 
atoms from the hot to the cold component (i.e., cool the 
atoms). The value of x^ is the largest ever reported, ex- 
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ceeding that obtained for EIT-based SWM by 10 7 
and in three-photon absorption in a zinc blend semicon- 
ductor by 10 24 [13. 
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For the same detuning (A = — 31"), we find Xx^yy — 
2.0(±2.3) x 10 -10 (m/V) 2 , where the large experimental 
uncertainty arises from our inability to accurately mea- 
sure /3 at the smallest I p . This x^ 3 ^ response is ~ 100 
times smaller than that observed in EIT-based systems 
[f|, and leads to four- wave mixing (FWM) in which the 
pump beams scatter off the weak atomic density grating 
composed solely of atoms in the hot component [35| . 

Beyond the region in which the Taylor series expansion 



given in Eq. [5] is valid 

,(5) 



, I c < I p <C I* t ), f c 1 and 
the x^ 5 - 1 response saturates. Nevertheless, /3 continues 
to increase linearly with I p , in contrast to x processes 
that increase only sub-linearly beyond saturation. 

We compare our observed NLO response to previously- 
reported nonlinearities based on x^ 3 ' processes by consid- 
ering the achievable NLO phase shift for a fixed I p . For 
example, Lo et al. Q observe a 0.25 rad phase shift for an 
intensity of 230 /jW/cm 2 . We find that /3L = 0.25 rad for 
I p = 560 /iW/cm 2 and A/T = —3, although our system 
does not require any auxiliary, strong coupling beams (as 
in the EIT-based setups). Also, j3 continues increasing 
quadratically with I p for our x^ process (rather than 
linearly as in x^ processes), which further aides us in 
achieving large phase shifts at low-light-levels. Our sys- 
tem provides the additional advantage that we can work 
far from resonance and thereby combine large nonlineari- 
ties with high transparency (similar to recent predictions 
for Rydberg-enhanced EIT |5|). We quantify the tradeoff 
between absorptive loss and NLO phase shift using the 
cross-phase-modulation figure of merit £ = /3/a& (i.e., 
the ratio of the NLO phase shift to power loss, see Fig. 
WX)- While Lo et al. Q observe a maximum value of 
C = 0.35 for an intensity of 230 /iW/cm 2 , we exceed this 
value for I p > 300 ^W/cm 2 and obtain a maximum value 
of £ = 26 for I p = 5.6 mW/cm 2 . Our system's ability to 
realize large £ arises directly from the independence of I c 
on detuning. 

For sufficiently large values of f3L, the back-action of 
the amplified probe fields influences the coupled light- 
matter dynamics, resulting in a strongly-coupled sys- 
tem with long-range atom-atom interactions. In this 
regime, additional cooling occurs in the x-z-plane due 
to the lattice formed by the amplified probe and nearly- 
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26j. We also observe 
c/10 5 (correspond- 



copropagating pump beams 
small group velocities of light v g 
ing to a slow light delay of td = 4 /is, see Fig. [3^,) by 
weakly modulating I s and recording the phase lag of the 
transmitted beam (relative to the case where I p = 0). 
The medium therefore acts like a high-finesse cavity and 
lengthens the photon lifetime in the gas such that the 
time scales of the atomic and optical dynamics become 
comparable (contrary to the assumptions leading to Eq. 
[2]). As a result of this back-action at large (3L, the 
atomic coherence time (defined as the 1/e decay time) 
increases to over 100 /is, which is over 50 times larger 
than that measured for f3L <C 1, where thermal atomic 



motion causes grating washout (see Fig. [3Jd). In addi- 
tion, we observe a collective instability for j3L > 1 in 
which the probe fields and density grating are sponta- 
neously self-generated via the NLO interaction. We will 
describe these results in a future publication. 
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FIG. 3. a) Dependence of the group velocity and slow light 
delay on the NLO coupling strength, b) Decay of the nonlin- 
ear coupling coefficient as a function of time after the signal 
beam is turned off at t = for A = — 5F We measure the 
long (short) decay time of 100 (1.8) i-ts using I p = 1.25 (0.45) 
mW/cm 2 , which corresponds to /3L= 0.84 (0.11). The inset 
shows a detailed view at small times. 



Our results are not limited to the particular experimen- 
tal geometry. We observe qualitatively similar results for 
other pump field polarizations (e.g., Iinf91in and <j + -<j + ), 
which indicates that the details of the configuration are 
not essential to the nonlinearity. On the other hand, we 
find that j3 ~ 100 x smaller for A > 0, where heating 
and limited bunching occur. Thus, the NLO process oc- 
curs for any setup that permits both atomic cooling and 
bunching, which is consistent with our physical interpre- 
tation of the nonlinearity. 
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